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ABSTRACT 

Vie discuss If and un<ter what conditions the coinblned emission from 
power>1aw sources can mimic the XRB spectrum In the range 3-50 keV measured 
with the A2 experiment on HEAO-1. We confirm that a good fit can be obtained, 
but the required spectral properties of component sources differ from those 
observed for local active galactic nuclei. Strong constraints are deduced for 
the lew luminosity extension and for the evolution of such local objects. We 
show that ai^ other class of sources significantly contributing to the X-ray 
background must be characterized by an energy spectral Index y < 0.4, the mean 
Index of the XRB (3-15 keV), and must exhibit steeper spectra at somewhat 
higher energies. 


Subject headings: galaxies: cluster of - radiation mechanisms 


X-rays: sources - X-rays: spectra 



I. INTRODUCTION 

Although almost every class of known extragalactic objects have been 
found to be X-ray emitters, It Is still unclear whether their Integrated 
emission can account for the full Intensity of the X-ray background (XRB). 

For the 2-10 keV band the most recent estimates yield a ~ 5% contribution from 
clusters of galaxies {McKee et al. 1980; Hintzen et al. 1980; Picclnottl et 
al. 1981) and a contribution » 20% for active galactic nuclei (Picclnottl et 
al. 1981). Moreover consideration of the constraints set by the deep survey 
counts (Glacconi et al. 1979), by the fluctuation level (Schwartz 1980) and by 
the optical counts at faint magnitudes (Kron 1980; Bahcall and Soneira 1980) 
has led to the conclusion that the contribution from quasars can not exceed 
50% (Cavallere et al. 1981). 

The 3-50 keV XRB spectrum is very accurately described by a 40 keV 
thermal bremsstrahlung njodel (Marshall et al. 1980). Known spactra of single 
discrete sources do not resemble that of the XRB. The X-ray enitt ig plasma 
In clusters of galaxies has in most cases a temperature < 10 keV (Mushotzky et 
al. 1978). Excluding BL Lac type objects, the 3-50 keV spectra of a sample of 
18 local active galaxies are power laws with energy Index 7 very close to 0.7 
(Boldt 1980) and no evidence of steepening of the spectrum up to at least 80 
keV has been found In the 8 objects for which higher energy measiirements have 
been obtained (Mushotzky 1980, private communication). As for quasars, our 
spectral knowledge Is limited to only two objects; although a thermal 
C'imponent compatible with the XRB Is permitted for 3C273 (Worrall et al. 

1977), radio quiet QSO 0241+622 requires y > 0.7 (Worrall et al. 1980), 

In order to account for the full XRB Intensity either some superposition 
of suitably evolving active galactic nuclei or a new astrophysical component 
Is apparently needed. While severe constraints are dictated by source counts 



and by the fluctuation level, complementary Information and additional 
constraints come from the observations of the XRB spectrum. A preliminary 
discussion, based on pre-HEAO-1 data, of the necessary conditions for the 
combined emission of non-thermal sources to mimic the observed spectral shape 
from 2 to 400 keV was carried out by Cavallere et al. (1979, Paper I). This 
analysis was preliminary In that It was necessary to combine the existing 
measurements from basically different experiments so that systematic errors 
dominated the results: this fact was made apparent by the large values 

O 

of X (> 2.^ per degree of freedom) which were obtained In any case (Including 
thermal bremsstrahlung). 

The A2 experiment on HEAO-1 has now provided a homogeneous set of high 
precision measurements, though covering a more limited energy range (3-50 
keV). We present here a new analysis of the constraints set by such data on 
spectral properties of power-law sources that In principle could make up the 
XRB. 

II. CLASSES OF MODELS 

The guidelines for the choice of the empirical classes of power-law 
spectral models to be explored are supplied by an Inspection of the observed 
XRB spectrum, which suggests the following strategy: 

(I) The bright resolved sources may (In fact, they likely do) dominate the 
XRB only at energies > 100 keV (Goldt 1981). Progressively fainter and more 
numerous sources are required to make up the XRB at progressively kwer 
energies; but then the component spectra should steepen or cut-off at 
characteristic energies E < 20 (1+z) keV to conform to the steepening of the 
XRB spectrum at « 20 keV. 

(II) As this svoepening appears now a gradual one with the present, sharply 
defined data, component spectra with a finite cut-off disbributlon are 



required to successfully describe the smooth transition across 20 keV (note 
the variance with the pre-HEAO-1 notion of a sharp break In the XRB spectrum 
that required sharp and fixed cut-offs of the component spectra). This 
requirement can be relaxed if E corresponds Instead to a spectral break that 
Is less drastic than a sharp cut-off. 

Only two main classes of non- thermal sources comply with the above 
criteria (cf. Paper 1): 


1) Sources whose power-law spectrum steepens at some energy Ej^; 

-Y. 


t(E) = t. 




E ^ E, 


E > E, 


( 1 ) 


2) Sources having a power-law spectrum with a sharp cut-off at 
an energy E^ : 
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In order to reproduce the XRB spectrum with Model 2 a broad distribution 
of cut-off energies must also be assumed. Two possibilities have been 
considered: 

2a) evolves with redshift from the local observed values 
Eco^ some 10* keV obeying 


( 3 ) 



2b) the sources have, at any z, a power law distribution of cut-offs: 
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In all cases a number/luminosity evolution «(l+z)“ was allowed for. 

The possibility that the XRB is saturated by an extrapolation at the 
faint end of the luminosity function of active galactic nuclei, without 
evolution, was also considered in Paper I. This possibility will not be 
discussed further here (except for a remark in Sect. 4) since it appears to be 
inconsistent with the flat faint end of that luminosity function as recently 
determined by Piccinotti et al. (1981). 

As in Paper I, we allow for a gaussian distribution of spectral indices 
with mean values y, y^, and dispersions o, and Og respectively. We 

have neglected, however, the dispersion of (model 1) or of E(.q (model 2a); 
this allows us to obtain (see Appendix), at least in the two 
cases a=l and n*0 (n being the density parameter), explicit analytic 
expressions for the combined emission 1(E) of sources and therefore to keep 
within reasonable limits the computer time needed for the subsequent analysis. 

III. RESULTS 

2 

As in Paper I we have applied the minimum x technique (see e.g. Avni 
1976; Cash 1976) to determine the best fit values and the allowed ranges for 
the parameters. Minimizations have been carried out by the MINUITS CERN 
subroutine package, written by James and Roos. 



The present ana1>sls confirms the result of Paper I, that a suitable 
superposition of power law sources can represent the X!^B spectrum as 
accurately as a thermal b^emsstrahlung distribution. This Is Illustrated In 
Figs. 1 and 2 , where the counts expected by convolving the spectrum predicted 
by Model 2a with the response functions of two detectors are compared with the 
observed ones. Table 1 gives, for all available data sets, the minimum 

9 0 2 2 

values of x » the mmibers of degrees of freedom v and Xy*x„,^n/v assoc- 
iated with both power law and thermal bremsstrahlung models. Since the spread 
In the results from different detector layers exceeds that expected from 
counting statistics alone. In computing these values a 2 % contribution of 
systematic errors to the overall uncertainty has been allowed for. We have 

not distinguished between power law models 1, 2a and 2b since the 
2 

corresponding Xj^^^ differ significantly. While It Is possible to fit 

the HEAO-1 data with these models, the fits are much more tightly restricted 
than they were by the composite data used for Paper I. It Is thus more 
significant that fits can be obtained. 

For the low energy slope and Its dispersion the results of fitting the 
HEAO-1 data differ from the best fit results In Paper I, although they He In 
the large possible ranges of parameters found there. No sign has been found 
In the HEAO-1 data of the slight flattening of the XRB spectrum from 3 to 15 
keV which Is to be expected In the presence of a large dispersion In the 
spectral Indices of contributing sources (the slope of 1(E) Is given 
by -Y + j o^ln E/2/J. cfr. eqs. (A1)-(A4)) and which seemed to be born out by 
pre-HEAO-1 data (cfr. Fig. 1 In Paper I). The HEAO-1 data require that the 
dominant contributors to the XRB have a rather narrow range of spectral 
Indices centered around y*0.4*0.5. This point Is specified In Fig. 3 which 
shows the 90% confidence contours for y (or Yj) and o (or o^) following from 



both HED1 and HED3 combined layers data (these contours are essentially 
Identical for all models considered here). As expected, the constraints on 
the slope of Model 1 are not so sharply defined. Our spectral fits 
Indicate 1.2 <y 2<1*4 with O2<0.35. 

The other Important parameter for comparison with direct observations is 
the energy E at which the spectral break or the cut-off occurs. For model 1 
we find a best fit value E|j » 50-55 keV (the exact value being slightly 
different when data from different detectors are used) with an allowed range 
(90% confidence) 40-60 keV. As for model 2a, the limited energy range covered 
by HEAO-1 data does not allow a precise definition of E^.^. The formal best 
fit value Is E^q <* 500 keV with a 90% confidence lower limit > 300 keV; 
the upper limit Is undetermined. The value of n Is strongly correlated with 
the value of E^.^ (It Increases with E^q following the approximate relationship 
given by Eq. (6) of Paper I): we have n*1.4 If E^.^ » 500 keV, n*l.l for 
Epo “300 keV, and n“2.25 for E^^»1500 keV. 

The quality of the fit turns out to be rather Insensitive to the assumed 
value of o (I.e. to the rate of evolution) as far as models 1 and 2b are 
concerned; In both cases however values of o>3 are favoured, lower values 
producing an Insufficient bending down of 1(E) because of the smearing out 
effect of the redshift distribution. For the same reason we expect, 
especially In the case of mdel 1, that a stronger evolution be associated 
with spectra having a round knee Instead of a cut-off or a sharp break. On 
the other hand, for model 2a the rate of evolution plays a direct role In the 
definition of the effective distribution of cut-off energies, see Eq. (3); 
only with 3.5 <a<5 can a good fit be obtained. 

IV. DISCUSSION 

If the XRB spectrum Is to be Interpreted In terms of emission from 



unresolved sources, Its remarkable smoothness In the 3-50 keV range seems to 
demand that most of the observed Intensity c(nnes from a single, fairly 
homogeneous population with a mean spectral Index Yj (or y) definitely flatter 
than that of local AGN's, and with a cut-off or a break occurring at 
relatively low energies, again at variance with the (however limited) 
observations of local AGN's. These two facts pose serious difficulties to all 
attempts at relating the main contributors to the XRB with local AGN's. For 
example, although model 2a Incorporates the large values of suggested by 
the presently available observations. It Is constrained by the observed local 
luminosity function for AGN's and requires an evolution In nunk>er/lum1nosity 
associated with a downward evolution of E^. Furthermore, to recover an 
effective distribution function of y'$ sharply peaked around T»0.4, the 
spectral Index of the XRB from 3 to 15 keV, the spectral Indices of the 
sources must also evolve (I.e. they must, on the average, decrease with z 
Increasing) and such evolution must be much faster than that of E^. For 
example, an exponential evolution of y (Y^'^^expC-qll-t/tQ)] with y^ • 0.7 and 
q <xl.3*1.4 coupled with an exponential luminosity evolution 
(L-L^expCOd- '/tg)] with Q » 4 seems to be needed to obtain a good fit of 
the < 15 keV XRB spectrum; but an exponential luminosity evolution also 
results In an effective distribution of E^ which Is not broad enough to avoid 
a steep decline of the computed XRB spectnnn at E > 15 keV, for any value 
of n (cf. eq. (3)). 

In the case of models 1 and 2b, continuity with local AGN's Is not 
required (and forcing such a continuity would anywqy bring In the previously 
noted difficulties). On the other hand. If the main contributors to the XRB 
are unrelated to AGN's the possible contribution to the XRB of the AGN's 
themselves must be < 30%, because a larger contribution from objects 



with y" 0.7 would prejudice any fit. This limit Is the most stringent to date; 
It Implies both that the luminosity function of AGN's with y« 0.7 must converge 
not far below U2-10) keV « 10^^ erq/s, and that essentially no evolution of 
these objects Is allowed, In X-rays at least. The only possibility for a 
downward extension of the AGN luminosity function would In^ly * T abruptly 
switching to a value » C.A at L » 10^^ erg/s. 

Indeed, the spectral properties of any population of power-law sources 
significantly contributing to the XRB Intensity are severely constrained. 

Since the allowed dispersion of spectral Indices Is already almost 

completely saturated by AGN's, any other population must be y (or y^) < 0.4, 
the mean spectral Index of the XRB In the range 3-15 keV. This Implies, for 
example, that QSO's either have the same spectral Index as the yet unknown 
population which saturates the XRB Intensity, or have a steep spectrum and 
make only a small contribution to the XRB above 3 keV. 

Similar considerations would constrair the evolutionary properties of 
sources (such as BL Lac's) having a very flot (y»0) spectral Index well above 
3 keV. On the other hand, the spectra of BL Lac type objects tend to be 
relatively steep at < 3 keV (Holt 1980) Indicating that their contribution to 
the XRB could be substantial only for the soft part of the spectrum that Is 
still not well defined. 
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Aop»ndlx 


Th« relevant formulaa for computini', tht eombintd •missions of 
power law sources ere given in the Appendix of Paper I to which 
the reader is also referred for nomenclature. Neglecting the di> 
spersionu of the break and cut'off energies (i.e. setting 
*^^b(c)^ * ^^^b(c) " ^'bCc)^ replacing (Its) with its weighted 

average expressions 0.5 o^lndts) (and in the 

analogous ones with o| and o’) we obtain the results listed be* 
lew for the various models. We have checked that the above men- 
tioned approximation introduces an error smaller than 1% for the 
interesting values of o's (o<0.5). 

For model 1 the integr«.^ed emission reads: 
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In actual calculations wso.S (corresponding to 0*1; sac aq# (5) 
of Paper I) and have bean adopter. The normalisation anar 

gy 2/$ kaV is fixed by our assumption of a 2-10 kaV luminosity 
function independent of the spectral index. 

For model 2a -c’.!* result is: 
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Finally for model 2b we have: 
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6 is the same as for the model 2a 
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Table 1 


Power law Thermal brem^strahlung 


Detector 

X ^ • 

V 


Xmin 

V 


HED 1 

Layer 1 

22 

17 

1.3 

24 

21 

1.1 

Layer 2 

23 

12 

1.9 

30 

15 

2.0 

Combined 

Layers 

33 

17 

1.9 

23 

21 

1.1 

HED 3 

Layer 1 

24 

17 

1.4 

30 

21 

1.9 

Layer 2 

11 

12 

0.90 

12 

15 

0.8 

Combined 

Layers 

14 

17 

0.85 

20 

21 

1.0 


MED 1 

19 19 


Layer 1 


1.0 





Figure Captions 


Fig. 1 - Comparison of the counts from layer 1 of HED 1 with tho 
se predicted by Model la with a=4.0, Y=0.465, 0=0.065, 

E =400 keV, n=1.35. 

CO 

Fig. 2 - Comparison of the counts from layer 1 of HED 3 with tho 
se predicted by Model 2a with the same values of the pa 
rameters as for Fig. 1 except for n=1.25. 

Fig. 3 - Constant y* contours at ■♦‘4.6 in the y-a plane 

for HED 1 (continuous line) and HED3(dashed line) comb^ 
ned layers data. If we do not care about the values of 
all other parameters, these contours bound the 90% confJ_ 
dence region for "y and O. The best fit values of the pa- 
rameters corresponding to HED 1 combined layers (O), HED 3 
combined layers (®) and MED layer 1 (0) are also shown. 
Model 2a. 
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